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Direct  Experimental  Evidence  for  an 
Atmospheric  Gravity  Wave  Cascade 

1.  INTRODUCTION 

The  saturated  cascade  (S-C)  theory  of  atmospheric  gravity  waves  (GW’s)  in  Dewan  (1994, 
1997)  utilized  both  “saturation”  and  “cascade”  concepts  to  predict  the  power  spectral  densities 
(PSD’s)  of  vertical  and  horizontal  wave  numbers,  and  temporal  frequencies  of  temperature, 
density,  and  wind  fluctuations  due  to  GW’s.  “  Saturation”  explained  the  quasi-universality  of  the 
vertical  wave  number  PSD’s  of  horizontal  wind,  temperature,  and  density  while  “cascade” 
explained  the  very  much  larger  variability  of  all  the  other  PSD’s.  The  combination  of  both 
saturation  and  cascade  concepts  also  resulted  in  an  explanation  of  the  observed  wave  period- 
wavelengths  seen  in  the  data  of  Gardner  and  Voelz  (1987),  Reid  (1986),  and  Manson  (1990).  At 
present,  no  direct  test  of  the  cascade  part  of  the  S-C  theory  has  been  performed.  This  would 
consist  of  the  simultaneous  measurement  of  both  the  average  turbulent  dissipation  rate,  s  ,  and 
certain  of  the  GW  PSD’s  strongly  dependent  on  s  (see  Dewan,  1997).  For  this  reason  it  is 
desirable  to  find  other  evidence  for  the  cascade  in  the  existing  experimental  literature.  Therefore, 
in  this  report  the  horizontal  wavenumber  PSD’s  of  wind  fluctuations  reported  in  Nastrom  et  al 
(1987)  and  Bacmeister  et  al  (1996)  are  examined  and  the  discontinuous  slopes  of  these  spectra 
are  explained  in  a  manner  which  is  in  quantitative  agreement  (so  far  as  estimation  permits 
without  further  data)  with  the  S-C  model  and  the  cascade  property  in  particular. 

The  cascade  described  in  Dewan  (1997)  starts  with  a  “  source”  at  the  large  scale,  low 
wavenumber  denoted  by  kz*.  The  growth  of  the  waves,  as  they  ascend,  provides  the  energy  per 

unit  mass  at  this  scale.  The  cascade  goes  from  this  large  scale  to  smaller  and  smaller  scale  (and 
period)  waves.  As  the  wave  field  ascends,  kz*  decreases  as  explained  in  Smith  et  al  (1987). 
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2.  EXPERIMENTAL  EVIDENCE 

Figure  1,  from  Nastrom  et  al  (1987)  shows  horizontal  wavenumber  spectra  from  horizontal 
velocity  profiles  based  on  the  well  known  GASP  program  measurements  made  from  a 
commercial  aircraft.  These  data  were  selected  to  ascertain  terrain  effects  and  do  not  reflect  an 
average  on  all  data  sets.  Our  Figure  1  (their  Figure  6)  was  further  restricted  to  the  “high  wind 
speed  cases”  over  mountains  (see  Nastrom  et  al  (1987))  and  was  further  segmented  to  reflect 
tropospheric  and  stratospheric  regions.  As  these  authors  pointed  out  there  are  transition  scales 
where  the  PSD  slopes  change  from  -5/3  to  -3  as  one  progresses  to  higher  wave  numbers.  These 
occurred  at  23  km  in  the  troposphere  and  15  km  in  the  stratosphere.  In  contrast,  Figure  2,  taken 
from  Bacmeister  et  al  (1996)  made  similar  measurements  from  an  ER-2  (also  called  a  U-2) 
aircraft  at  20  km  altitude.  (This  is  their  Figure  14b).  They  also  observed  a  change  from  -5/3  to 
-3  with  the  steeper  slope  associated  with  scales  smaller  than  3  km.  (See  their  abstract.)  The 
three  observations  of  the  scale  of  the  breaks,  which  occurred  at  wavelengths  of  25  km  in  the 
troposphere,  15  km  in  the  stratosphere,  and  3  km  much  higher  in  the  stratosphere  suggest  that,  in 
accordance  with  Dewan,  1997,  a  cascade  (showing  a  -5/3  slope)  takes  place  from  higher  scales  to 
lower  scales  as  the  gravity  waves  propagate  upwards  in  altitude.  The  main  sources  are  assumed 
to  be  located  below  the  regions  of  observation.  This  phenomenon  was  predicted  and  noted  in 
Dewan  (1997)  and  is  based  on  the  “cascade  set  up  time”  hypothesis  of  W.  Hocking,  (private 
communication,  1996). 

In  other  words  the  spectral  break  occurs  because  as  the  cascade  is  “  set  up” ,  that  is,  extends  to 
sequentially  smaller  scales,  it  takes  an  extended  time  before  it  reaches  the  smallest  wave  scales. 
Note  that  previous  theories  of  the  break  label  the  break  wavenumber  as  k*.  It  is  important  to 

avoid  this  confusion.  Note  also  that  k*  decreases  with  altitude  in  contradiction  to  the  observation 

that  the  k  of  the  break  increases  with  altitude,  thus  disproving  these  previous  theories. 

It  should  be  mentioned  that  the  majority  of  experimental  horizontal  wavenumber  spectra  of 
winds  do  not  contain  breaks  from  -5/3  to  -3  slopes  (see  Nastrom  and  Gage  (1985)  and  U-2 
stratospheric  observations  (20  km)  of  Crooks  et  al  (1968)).  Instead,  the  breaks  seem  to  occur 
only  when  there  is  an  obvious,  relatively  near,  source  of  gravity  waves  as  in  the  case  of  Nastrom 
et  al  (1987).  In  future  experiments  this  can  be  investigated  further.  Additionally  one  would  also 
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expect  that  the  presence  of  breaks  would  be  seen  more  often  at  high  altitudes  since,  for  example, 
doubling  the  altitude  would  quadruple  the  area  of  exposure  to  wave  sources  below. 

The  specific  goal  here  is  to  show  that  the  time  required  for  a  GW  cascade  field  to  ascend 
from  13  km,  i.e.  the  nominal  stratospheric  altitude  in  Nastrom  et  al  (1985)  (private 
communication  with  Nastrom,  1998),  to  the  higher  stratospheric  altitude  (20  km)  of  Bacmeister 
et  al  (1996)  is  about  equal  to  the  estimated  time  needed  for  the  cascade  to  progress  from  the  15 
km  lower  stratospheric  scale  to  the  3  km  upper  stratospheric  scale.  This,  as  will  be  shown, 
therefore  quantitatively  explains  the  scale  of  discontinuous  change  in  PSD  slopes  and  its  altitude 
dependence  as  shown  here  in  Figures  1  and  2. 


Zonal  Wind  Spectra 
High  Wind  Speed  Cases 


Figure  1.  Horizontal  wavenumber  spectra  of  zonal  winds  (high  wind  speed  case  over  mountains) 
from  tropospheric  and  stratospheric  measurements.  Note  the  scales  of  the  breaks  where  the 
slopes  change  from  -5/3  to  -3.  This  figuure  is  from  Nastrom  et  al  (1987)  and  it  is  their  Figure  6. 
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Figure  2.  Horizontal  wavenumber  spectrum  of  horizontal  winds  at  20  km  altitude  from 
Bacmeister  et  al  (1996).  It  is  their  Figure  14b.  There  is  a  break  from  a  slope  near  -5/3  to  one 
significantly  larger  at  around  3  km. 


3.  QUANTITATIVE  TREATMENT  USING  THE  S-C  MODEL 

The  quantitative  model  to  be  described  requires  the  knowledge  of  a  number  of  parameters. 
Since  they  are  not  presently  available  it  will  be  necessary  to  use  nominal  values  based  on  past 
observations.  They  could  be  obtained  in  future  experiments.  These  are  listed  in  Table  IB  (based 
on  Dewan  (1997)).  In  other  words,  only  these  already  published  values,  rather  than  values 
“tuned”  to  Figures  1  and  2,  are  used.  Our  present  purpose  is  merely  to  see  if  the  model  gives 
results  that  are  within  reason.  Table  1A  also  gives  the  required  equations  from  Dewan  (1997). 


Table  1. 


A.  Equations  (from  Dewan  (1997),  Table  1).  (See  B  below  for  constants). 

Horizontal  velocity  horizontal  wavenumber 

spectrum 

Vertical  velocity  horizontal  wavenumber 

spectrum 

<oo-g  ’fNVa|)<"  P-)K“ 

Horizontal  velocity  frequency  spectrum 

(  x  a  •  (s  aj/a,)  x  -2 

Vvx(©)_  2  l'(2n)a> 

Characteristic  vertical  wavenumber 

Saturated  cascade  relation 

K 

(saj/a,) 

Wavenumber  bandwidth  (from  Eq.  (21), 

Dewan  (1997)) 

Ak,-A 

4n 

Wavelength-period  relation 

x3/2  =  f=2s_y'2  ^ 

vsaj/a,^ 
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B.  Parametric  values  (Dewan  (1997),  Table  2) 


Stratospheric  dissipation  rate  (volume  and  time  s  = 
average) 


Mesospheric  dissipation  rate  (volume  and  time  g  = 


average) 


Brunt  frequency 
Brunt  period 


Characteristic  vertical  wavenumbers  based  on 
Smith  et  al  (1987) 


N  =  0.02  rad/s 


xB  =  =  3i4s 

B  N 


43°  latitude  (representative)  inertial  frequency  f  =  10-4  rad/ 


ropospneric  k  =zn  /  i  ion 

z* 

Stratospheric  k  =  2n/5  km 
z* 


The  cascading  GW  field  in  Dewan  (1994, 1997)  originates  with  the  largest  horizontal 
wavelengths  (Xx)  and  periods  (x)  and  cascades  downwards  to  the  small  wavelengths  and  periods. 
An  estimate  of  the  vertical  ascent  rate  of  this  GW  field  would  be  of  the  order  of  the  vertical 
group  velocity  ugz  of  the  inertial  waves  (see  Dewan,  1997).  This  ugz  is  sufficiently  well 

approximated  by 


o„  =  —L 


where  Xz*  is  the  longest  S-C  wavelength  or  the  “  characteristic”  wavelength  (c.f.  Smith  et  al 

(1987)  and  Dewan  (1997))  and  xf  is  the  inertial  period.  Let  TOF  equal  the  “time  of  flight”  or 
ascent  time  for  the  GW  field  to  go  from  13  km  to  20  km.  Then 
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(2) 


TQP  _  (20  km  -  13  km) 

u„ 

»z 

At  43°  Lat.,  Tf  =  =  6.28  x  104s  approximately  from  Table  IB.  The  nominal  value  for  in 

the  stratosphere  from  Smith  et  al  (1987)  and  Table  IB  is  5  km.  Using  these  in  Equation  (1)  and 
Equation  (2)  gives  ogz  =  8  x  10-2  m/s  and  TOF  =  24.3  hrs  respectively. 

We  next  compare  this  to  the  “  cascade  duration” ,  CD,  needed  to  propagate  down  in  scale 
from  A.x  =15  km  (seen  at  13  km  altitude)  to  Ax  =  3  km  (seen  at  20  km  altitude).  The  following 
calculations  depend  crucially  on  information  in  Dewan  (1997)  and  the  reader  must  refer  to  Table 
1 A  for  the  equations  used  below.  Thus  the  cascade  duration  is  given  by 

J 

CD  =  £  nxj 

i=0  (3) 

where  tj  is  the  period  of  the  ith  wave  mode,  n  is  the  number  of  periods  needed  for  one  mode  (e.g. 
the  (i  -  l)th  mode)  to  decay  and  to  simultaneously  “feed”  the  next  smaller  mode  (e.g.  the  ith 
one),  and  J  is  the  total  number  of  “cascade  steps”  (in  the  model)  undergone  during  the  cascade 
duration  under  investigation  (c.f.  Dewan  (1997)).  In  other  words  J  is  the  number  of  modes 
involved  in  going  from  the  initial  (larger)  scale  to  the  final  (smaller)  scale  as  the  wave  field 
ascends. 

Table  1 A  contains  the  wavelength-period  relation 

/  \2/3 

Ti  =  (AXx)  (4) 

where  A,x  refers  to  horizontal  wavelength  and  where 

(5) 
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where,  again,  s  is  the  time  and  volume  average  turbulent  dissipation  rate,  and  the  ratio  (aj/a2)  is 
a  universal  constant  in  the  model.  (See  Table  1.)  From  Table  1,  a2/aj  =  2.5  x  103.  Also,  from 
Table  1  an  appropriate  nominal  average  value  for  stratospheric  s  is  10-5  m2/s3.  (It  would  be 
highly  desirable  that  average  s  be  measured  simultaneously  with  the  PSD’s  in  question  in  a 
future  experiment.)  These  nominal  values  give  one  A  =  15.9  which  should  be  regarded  as  a 
“  rough”  estimate  to  be  tolerated  on  the  basis  that  it  is  the  best  available  for  the  present.  Thus 
one  can  estimate  Xj  for  any  given  XXj  from  Equation  (4). 

The  next  step  in  our  estimation  of  CD  from  Equation  (3)  is  to  obtain  J  as  well  as  the  values  of 
the  intermediate  X\  s  (and  xj  s)  of  each  step  in  the  model  cascade.  For  this  we  start  with  the  mode 
bandwidths  given  in  Table  1 A  by 


A  kz  -  ^  (6) 

where  kz  is  the  vertical  wave  number,  2it/Xz,  and  where  Akz  is  the  bandwidth  of  the  mode  in 
question.  However,  Equations  (3)  and  (4)  require  X,Xj  and  AA,Xi  (i.e.  the  bandwidth  in  terms  of  the 

horizontal  wavelength). 

For  this  purpose  we  utilize 


K 


^a,  N3^ 


1/3 


V  a2s  J 


.  1/3  ^  1/3  ,  1/3 

kx  =  G  kx 


(7) 


where  kx  is  2nlXx  and  where  Equation  (7)  is  given  in  Table  1  above.  We  use  Equation  (7)  in 
order  to  relate  A-Xj+1  to  A,Xj  (i.e.  the  wavelengths  of  contiguous  modes  separated  by  bandwidth 

AlXj).  We  start  by  defining 


xi+l 


=  A,v.  +  AXV 


(8) 


8 


In  other  words,  AAXj  is  the  cascade  step  size.  This  relation  can  be  deduced  from  Equations  (6) 
and  (7).  Differentiating  Equation  (7)  and  using  Equation  (6)  yields 


Akz=  3 


rl/3  9  r 

—  K  Ak, 


4n 


(9) 


Solving  the  second  equality  in  Equation  (9)  for  Akx  and  using  Equation  (7)  to  eliminate  kz  from 
Eq.  (19)  gives 


Akx 


3k, 

4n 


To  obtain  AAX,  one  differentiates  kx  =  2n/Xx  obtaining 


(10) 


Akx  =  ??.(-■)  AXx 
x  2 

Hence,  from  Equations  (10)  and  (11) 

Ak  =  ^x  _  (271)  3 
X  '  4"  "  \(4n)  ' 

Therefore,  solving  Equation  (1 1)  for  AAX  and  using  Equation  (12) 


(11) 


(12) 


AX 


x 


(13) 


Using  Equation  (13)  in  Equation  (8)  gives 


+  AX 

xi 


(14) 


or 
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where  i  =  0  refers  to  the  scale  at  the  start  of  the  CD  (that  is  to  say,  the  largest  break  point  scale  or 


A,max)  and  ^Xj  is  the  scale  under  consideration.  The  smallest  scale  would  be  A,mjn  or  Aj  the  scale 

reached  at  the  end  of  the  CD.  Thus,  in  our  example  A,mjn  or  Aj  =  3  km  and  A,max  or  A,o  =  15  km. 
Using  Equation  (15),  then,  we  obtain 


gJ  _  ^min 


max 


(16) 


Hence  one  can  obtain  the  value  of  J  from  Equation  (16)  as 


t  _  i*1  (^'mir/^'max) 

ln(B) 


(17) 


One  can  now  evaluate  Equation  (3)  for  the  CD  by  using  x\  from  Equation  (4) 


J  J  /  X 

CD  =  X  ntj  =  nE  (A Ax.) 
i=0  i=0 


2/ 


(18) 


and  using  Equation  (15)  for  A.Xj 


CD  =  a  X,  (AB  'iJ’-nA* 

i=0  i=0 


(19) 


or 


CD  =  „A“  (O28  2  (b2'T 

i=0 


(20) 


The  sum  in  Equation  (20)  is  readily  evaluated  by  using  the  general  formula  for  the  sum  of  a 
geometric  series 


M 

I  ar‘  = 


i=0 


a 


lr-1 


(21) 


(see  Burington  (1954)  p.  3).  This  gives 


Inserting  this,  the  value  of  A  from  Equation  (5),  and  the  value  of  B  from  Equation  (14)  we  finally 
obtain  from  Equation  (20) 


CD  =  n 


aj  270 

a 2  8  J 


1/3 


2/3 


'"max 


(23) 


To  evaluate  this  we  use  Xmax  =  15  km,  Xmjn  =  3  km  and,  from  Table  1,  n  =  4.  This  yields  B  = 
0.813  from  Equation  (14)  and  J  =  7.77  from  Equation  (17).  Thus  Equation  (23)  yields,  using  the 

-  m2 

above  numbers  and  the  earlier  quoted  values  of  a2/ai  =  2.5  x  103  and  s  =  10  — , 

s2 

CD  =  22  hrs  .  This  is  to  be  compared  to  TOF  of  Equation  (2)  since,  from  theory,  CD  =  TOF  . 

As  will  be  recalled  TOF  =  24.3  hrs  which,  in  view  of  the  approximate  nature  of  the  calculation,  is 
too  good  an  agreement  to  be  expected  and  must  therefore  be  considered  partially  fortuitous.  If 
A.min  =  2  km  had  been  used,  CD  =  24.5  hrs. 
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Figure  1  suggests  an  additional  check  on  our  model.  Let  us  compare  TOF  and  CD  for  the 
passage  from  tropospheric  height  (10.5  km  Nastrom,  private  communication)  to  stratospheric 
height,  13  km.  The  Lmax  =  23  km  and  Xmin  =  15  km.  Using  the  same  parametric  values  as  before, 
s  =  10'5  m2/s3,  n  =  4,  a2/a,  =  2.5  x  103,  one  finds  CD  is  about  1 1  hrs.  As  seen  in  Table  1  the 
values  of  Lz<  are  nominally  1  km  in  the  troposphere  and  5  km  in  the  stratosphere  (we  don’t  have 

tropopause  height  information  or  actual  Lz>  information)  yielding  vg  =  1.59  x  10'2  m/s  and 

7.9  x  10'2  m/s  respectively.  The  resulting  TOF’s  are  44  hrs  and  9  hrs  respectively  thus  bracketing 
the  CD  of  1 1  hrs.  The  agreement  is  therefore  within  reason. 

It  should  be  mentioned  that  Gao  and  Meriwether  (1998)  also  measured  horizontal 
wavenumber  spectra  by  aircraft.  While  Nastrom  et  al  (1985)  (1987)  measured  only  horizontal 
velocity,  Bacmeister  et  al  (1996)  and  Gao  and  Meriwether  (1998)  also  measured  the  vertical 
velocity.  In  Bacmeister  et  al  (1996)  the  break  in  both  horizontal  velocity  and  vertical  velocity 
occur  at  the  same  scale.  This  could  be  explained  by  saying  that  this  is  the  scale  reached  by  the 
cascade.  In  Gao  and  Meriwether  (1998),  in  contrast,  there  is  no  break  at  all  for  the  horizontal 
winds  but  there  is  one  for  the  vertical  winds.  The  present  model  does  not  explain  this,  thus  one 
has  a  new  unanswered  question  requiring  further  study. 

4.  FUTURE  MEASUREMENTS 

To  test  the  model  suggested  here,  the  following  measurements  would  be  useful.  The  PSD’s 
should  be  measured  at  high  and  low  altitudes  at  nearby  times  and  places  and  under  similar 
conditions.  The  value  of  s ,  i.e.  the  time  and  volume  averaged  turbulence  dissipation  rate  at  the 
altitude  of  the  aircraft  should  be  measured.  It  would  also  be  useful  to  measure  XZf  (c.f.  Smith  et 

al.  (1987)).  The  presence  of  gravity  wave  sources  such  as  jet  streams,  thunderstorms,  and  wind 
over  mountains  should  also  be  documented.  In  this  way  one  could  verify  the  predicted 
association  between  spectral  break  and  GW  source.  Of  course  the  direct  tests  of  the  saturated 
cascade  hypothesis  suggested  in  Dewan  (1997)  should  also  be  very  useful  if  not  of  primary 
value. 
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5.  CONCLUSION 

Experimental  data  published  by  Nastrom  et  al  (1987)  and  Bacmeister  et  al  (1996)  give  some 
quantitative  support  to  the  cascade  hypothesis  as  contained  in  the  S-C  model  in  Dewan  (1997). 
Further  tests  of  this  type  would  entail  simultaneous  measurements  of  PSD’s  and  s  ’s  over  similar 
locations  and  conditions. 
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